Topological superfluid phases of an atomic Fermi gas with in- and out-of-plane Zeeman fields and 
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We analyze the effects of in- and out-of-plane Zeeman fields on the BCS-BEC evolution of a Fermi gas 
with equal Rashba-Dresselhaus (ERD) spin-orbit coupling (SOC). We show that the ground-state of the system 
involves novel gapless superfluid phases that can be distinguished with respect to the topology of the momentum- 
space regions with zero excitation energy. Depending on the combination of Zeeman fields and SOC, the zeros 
may correspond to one or two doubly-degenerate spheres, two or four spheres, two or four concave spheroids, or 
one or two doubly-degenerate circles. Such changes in the topology signal a quantum phase transition between 
distinct superfluid phases, and leave their signatures on some thermodynamic quantities. 
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Introduction. While the topological matter and related phe- 
nomena have become a popular subject especially in the past 
few years within the condensed-matter community fl}[3), a 
new window has just recently been opened to the exotic world 
of topological phases in atomic and molecular physics com- 
munities, with arguably a greater premise on a route toward 
studying them under highly controllable atomic settings |@}J7). 
The current advances in atomic systems offer the possibility 
of engineering both Abelian or non-Abelian artificial gauge 
fields on demand, by coupling the internal states of atoms to 
their center-of-mass motion via Raman dressing of atomic hy- 
perfine states with laser fields. While there are a number of 
theoretical proposals for implementing atomic gases with var- 
ious SOC symmetries, several experimental groups have so far 
only achieved an Abelian ERD coupling, first with Bose [4, 5 1 
and then with Fermi [6, 7 | gases. Allured by the experimen- 
tal possibilities, there has been a growing interest in study- 
ing spin-orbit coupled atomic Bose and Fermi gases as func- 
tions of the tunable laser parameters including the strength 
and symmetry of the SOC field, s-wave scattering length and 
Zeeman fields. Both the uniform and trapped systems have 
been considered in all one, two and three dimensions, already 
revealing exotic single-, two-, few- and many-body properties 
at both zero and finite temperatures ll8 W2"Tll . 

In particular to the Fermi gases, depending on both the 
strength and the symmetry of the SOC field, it has been shown 
that the interplay between the SOC and out-of-plane Zee- 
man fields gives rise to a variety of topologically-nontrivial 
and thermodynamically-stable gapless superfluid phases lTT~2l — 
|2D . In confined geometries, some of these phases have 
been shown to feature Majorana zero-energy bound states, 
for which the associated quasiparticle operators are self- 
Hermitian, i.e. a Majorana quasiparticle is its own anti- 
quasiparticle ED . These quasiparticles play a key role 
in fault-tolerant quantum computation, and although they 
are predicted to be found in strongly-correlated systems in 
many fields of physics, including the fractional quantum 
Hall systems l22l . chiral two-dimensional p-wave super- 
conductors [23, 24|/superfluids 11251 1261 . three-dimensional 
topological insulator-superconductor heterostructures 07], 
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FIG. 1: (color online) The gapless superfluid phases can be classi- 
fied with respect to the topology of their momentum-space regions 
with zero excitation energy. The top (bottom) figures illustrate the 
possible phases on the BCS (BEC) side, where the zero-energy re- 
gions correspond to (a,e) two (one) doubly-degenerate spheres cen- 
tered at k = when a = 0, ho 7^ and/or hi 7^ 0; (b,f) four 
(two) spheres centered at k x = iMa and k y = k z — when 
a 7^ 0, ho =0 and hi 7^ 0; (c,g) four (two) concave spheroids 
when a / 0,/io 7^ and hi 7^ 0; and (d,h) two (one) doubly- 
degenerate circles centered at k x = when q / 0,/io ^ and 
hi — 0. Here, a, ho and hi are the strengths of the SOC, out-of- 
and in-plane Zeeman fields, respectively. 



one-dimensional nanowires l28l [29], spin-orbit coupled 
semiconductor-superconductor heterostructures |30], etc., 
their realization has proved to be quite difficult. 

Given that the cold atoms offer unprecedented control in 
comparison to the condensed-matter ones, there is a good 
chance of creating and observing Majorana quasiparticles 
with atomic systems in the near future. Motivated by these 
recent developments |6] |7J, here we analyze the effects of 
not only the out-of-plane but also the in-plane Zeeman fields 
on the phase diagram of a Fermi gas with ERD coupling. 
Since the inclusion of an in-plane Zeeman field complicates 
the many-body problem considerably even at the mean-field 
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level, this term has been completely ignored in the recent lit- 
erature Ifl2l - ET1 . even though it appears naturally in the active 
experiments j6H2|. By including it here, we find gapless su- 
perfluid phases that can be distinguished with respect to the 
topology of the momentum-space regions with zero excita- 
tion energy, e.g. one or two doubly-degenerate spheres, two 
or four spheres, two or four concave spheroids, or one or two 
doubly-degenerate circles (see Fig. [TJ. These phases are dis- 
tinct because, even though their order parameters all have the 
same s-wave symmetry, changes in the topology leave strong 
signatures on some of the thermodynamic quantities, signal- 
ing a topological quantum phase transition. 

Mean-Field Theory. These results are obtained by consid- 
ering the following fields: the ERD coupling —ak x a y , out-of- 
plane Zeeman field —ho<J z and in-plane Zeeman field —hj<j y , 
where {a, ho, hi} > are their strengths. These fields have 
very different origin in cold-atom systems, where the Raman 
coupling and detuning between the two laser beams driving 
the Raman transition from the two-photon resonance corre- 
spond to the out-of and in-plane Zeeman fields, respectively, 
and the ERD coupling is the Doppler shift an atom experi- 
ences as it moves in the two laser fields 0-0. 

The mean-field Hamiltonian for this system (in units of H = 
1 = k B ) can be written as, H = H + (1/2) J2\c ^t^Lipk, 
where H = ^ k £k + | A| 2 /<7 and the matrix 
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A k f "k i) u, -k,f> °-k,J denotes the fermionic 
operators collectively, £k = £k — A* w i tn £ k = fc 2 /(2M) and 
fi the chemical potential, A = g(ak,t a -k,4.) is the mean-field 
order parameter with g the strength of the contact interaction 
and (• • • ) the thermal average, and S k = iak x . We eliminate 
g in favor of the s-wave scattering length a s via the relation, 
1/g = -MV/{Aira s ) + Y, k l/(2e k ), where V is the volume. 
The thermodynamic potential for this Hamiltonian becomes, 
fl = H Q + (T/2) J2k.x In [1 - f{Ek,x)} , where T is the tem- 
perature, A = {1,2,3,4} labels the quasiparticle/quasihole 
excitation energies £k,i = — -E-k,3 and _Ek,2 = — ^-k,4 
as determined by the eigenvalues of the matrix given above, 
and f(x) — l/(e x / T + 1) is the Fermi function. The self- 
consistency equations are obtained by setting dfl/d\A\ = 
and N = —dQ,/dp,, i.e. 
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Similarly, we obtain the in- and out-of -plane number polar- 
izations via Pj o = —(l/N)dn/dhi o, leading to Pj o = 
[-l/(2N)]J2 k>x (9Ek,x/dh I ,o)f(E k ,x). We also checked 
the stability of our mean-field solutions for the uniform su- 
perfiuid phase using the curvature criterion, i.e. a nonuniform 



superfluid phase (e.g. a phase separation) is favored when 
<9 2 fV<9|A| 2 < o. 

Excitation Spectrum. In general, all four branches of the 
excitation spectrum may be different from each other, and 
hence both £k,A and their |A|,/x, hj and ho derivatives are 
not analytically tractable and need to be evaluated simulta- 
neously with the self-consistency Eqs. |2| and (|3j for all k- 
space points. This makes the current work numerically more 
demanding compared to the previous ones [ 12-20 1, for which 
the explicit forms of i?k,A can be found. In order to under- 
stand the topology of the resultant superfluid phases, next we 
analyze the analytically- tractable limits. 

(i) No SOC (a — 0) limit. In this limit, we obtain 
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where si = s 2 = P2 = P3 = +1 and s 3 = s 4 = pi = p 4 = 
— 1 and hx = \/h 2 + hj is the strength of the total Zeeman 
field, showing that E^ i = —E kt 3 can have zero-energy re- 
gions in k space. The doubly-degenerate zeros are determined 
by h\ — £k — |A| 2 , and when ft/r > |A| these conditions give 
two spheres of zeros for [i > Wh^, — |A| 2 and one sphere of 
zeros for /i < ^Jh\ — |A| 2 . Therefore, the transition from 
gapped superfluid to gapless superfluid occurs at Ht — |A|, 
and /j, = determines the transition from gapless superfluid 
with two doubly-degenerate spheres of zeros to the one with 
one doubly-degenerate sphere of zeros, i.e. it gives the crit- 
ical point for the disappearance of the doubly-degenerate in- 
ner sphere. These possibilities are schematically illustrated in 
Figs.[TJa) and[TJe), and the corresponding gapless superfluid 
phases in this case are known to be topologically trivial lfT2l — 
ED. 

(ii) No out-of -plane Zeeman field (ho = 0) limit. Since the 
SOC breaks the inversion symmetry in the k x direction, we 
obtain 



Eu. 
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showing that E^, 



-ELk.a can have zeros. The zeros are 



determined by fc 2 + k\ = 2M(p ± ^Jh\ - |A| 2 ), where 
k x = k x — Ma, k± — ^Jk 2 + Ac 2 and ju = fj, — Ma 2 /2, 

and when hi > \A\ these conditions give four spheres of 
zeros for /I > \Zhf- — |A| 2 and two spheres of zeros for 
/x < \/h 2 j — | A| 2 . Therefore, the transition from gapped su- 
perfluid to gapless superfluid occurs at hj = |A|, and ju = 
determines the transition from gapless superfluid with four 
spheres of zeros to the one with two spheres of zeros, i.e. it 
gives the critical point for the simultaneous disappearance of 
two inner (one each from £k,i and E'k.a) spheres. Note that 
since the SOC gauge field can be integrated out via a shift in 
the k x momentum and fi, the corresponding gapless superfluid 
phases have the same topology and number as case (i), and 
therefore they are also topologically trivial. However, in com- 
parison to case (i), since the SOC breaks the inversion symme- 
try in the k x direction, a^O removes most of the degeneracy 
of zeros except for one or four circles, as schematically illus- 
trated in Figs.[TJb) and[TJf). Note also that sufficiently large a 
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splits these spheres of zeros completely away from each other, 
leaving no degeneracy. 

(Hi) No in-plane Zeeman field (hi = 0) limit. In this limit, 
we obtain 

£k,A = s x ^& + h 2 Q + a 2 k 2 x + |A| 2 + 2s xPx A ]il (6) 

where A k = V l o(£k + + " 2fc ^k> showing that 
Ek i = — i?k,3 can nave zeros. The doubly-degenerate ze- 
ros are determined by k x = and h 2 Q — ^ = | A| 2 , and when 
ho > |A| these conditions give two circles of zeros for fi > 
\/h 2 — |A| 2 and one circle of zeros for fi < \/h 2 Q — |A| 2 . 
Therefore, the transition from gapped superfluid to gapless 
superfluid occurs at ho = |A|, and [i = determines the 
transition from gapless superfluid with two doubly-degenerate 
circles of zeros to the one with one doubly-degenerate circle 
of zeros, i.e. it gives the critical point for the disappearance 
of the doubly-degenerate inner circle. These possibilities are 
shown in Figs.[TJd) and[TJh), and in sharp contrast to the cases 
(i) and (ii), the corresponding gapless superfluid phases in this 
case are known to be topologically nontrivial lfT4l[T6l . 

In the light of these limits, it is easier to understand the 
topology of zeros in the most general case when all three 
gauge fields are nonzero. This generic case is schematically 
illustrated in Figs.[TJc) and[TJg), showing that ho 7^ breaks 
the remaining degeneracy of one or four circles of zeros that 
are still present in Figs. [TJb) and[TJf). Therefore, the topol- 
ogy of gapless regions correspond to four concave spheroids 
on the BCS side and two concave spheroids on the BEC side, 
and the transition between the two occurs around /i ~ 0, with 
the simultaneous disappearance of the two inner (one each 
from E^ i and £^,3) concave spheroids. For instance, Fig. [2] 
shows the full k dependence of the momentum distributions 
for these phases, where kp = (S^N/V) 1 ^ 3 is the Fermi mo- 
mentum. It is important to emphasize that the exotic super- 
fluid phases discussed above, in particular the ones illustrated 
in Figs.[TJb),[TJc),[TJf) and[TJg), as well as the numerous pos- 
sible ways of transitions in between are unique to this work. 
To explore the feasibility of observing these phases in cold- 
atom experiments, next we analyze the phase diagram of the 
system at T = 0. 

Ground-state Phase Diagrams. The generic phase dia- 
grams are shown in Fig. [3] We note that, since the SOC can be 
integrated out of the self-consistency Eqs. (|2]) and Q on the 
x axis when ho = 0, the critical hi or Pi for the transitions 
from a gapped superfluid to a gapless one with two or four 
spheres of zeros and from the latter to the normal phase are ex- 
actly the same as those obtained for a population-imbalanced 
Fermi gas without the SOC, i.e. for the transition from a 
gapped superfluid to a gapless one with one or two doubly- 
degenerate spheres of zeros and from the latter to the normal 
phase. Here, we locate the normal phase boundaries by find- 
ing I A I < lCP 4 e^ where ep = kp/(2M) is the Fermi energy, 
and the inaccessible regions are determined by Pf + Pq > 1. 
In addition, on the y axis when hi = 0, our results recover the 
recent works |[T4l [161 . 

Away from the x and y axes, we find that the gapless su- 
perfluid phases with two or four concave spheroids of zeros 




FIG. 2: (color online) The momentum distributions nk,cr are shown 
as a function of k x and k± — for the superfluid phase 

with (a) two and (b) four concave spheroids of zeros. Here, we set 
l/(kpa s ) = and a = 0.5kp/M, and choose Po = 0.5 and 
Pi = 0.25 in (a), and P Q = 0.25 and Pi = 0.5 in (b). 



occupy quite large regions in the phase diagrams as shown in 
Fig. [3] Depending on a and a s , note that the transition from 
a gapless superfluid phase with two or four spheres of zeros 
to the one with two or four concave spheroids of zeros may 
require a finite threshold for ho- For instance, when the SOC 
is weak so that the spheres of zeros contain doubly-degenerate 
circles as shown in Figs.[TJb) and[TJf), an arbitrarily small ho 
splits these degeneracies, immediately giving rise to the tran- 
sition. This is clearly seen in Fig.|3jc). On the contrary, when 
the SOC is sufficiently strong so that it splits the spheres of 
zeros completely away from each other leaving no degener- 
acy, the transition occurs at a finite ho, as shown in Figs.[5Ja) 
andgfb). 

Next we note an important difference between the grand- 
canonical (Zeeman fields) and canonical (number polariza- 
tions) phase diagrams on the BCS side of the BCS-BEC evo- 
lution, and consider the ho — limit for the sake of its sim- 
plicity. In a grand-canonical phase diagram, where hi is in- 
creased from 0, it can be shown that while |A| is unaffected 
by hi 7^ until hi = |Ao|, where |Ao| is the value of 
the field-free order parameter when hi = 0, it immediately 
vanishes for hi > |A |. Therefore, increasing hi leads to 
a first-order quantum phase transition between a gapped su- 
perfluid and normal phase, without any intermediate gapless 
superfluid phase in between. On the other hand, in the case 
of a canonical phase diagram, where Pj is increased from 
0, there exists a |A| 7^ solution with h c < hi < |An|, 
where h c is determined by the maximum possible number po- 
larization Pi = P c . Therefore, in sharp contrast to the grand- 
canonical case, it is possible to obtain a gapless superfluid 
phase with increasing Pj. For this reason, the canonical phase 
diagram suits better for presenting the possible gapless super- 
fluid phases at unitarity as shown in Fig. |3jc). However, on 
the BEC side of the BCS-BEC evolution, such a complication 
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FIG. 3: (color online) The ground-state phase diagram is shown as 
a function of out-of- and in-plane (a,c) number polarizations and (b) 
Zeeman fields, where l/(&Fa s ) = 1 in (a, b) and in (c). Here, 
a — 0.5/cf/M in all figures, and the labels correspond to an In- 
accessible region (I), normal (N) phase, gapped superfluid (SF), and 
gapless superfluid with zeros of one or two doubly-degenerate circles 
(1DDC, 2DDC), two or four spheres (2S, 4S), or two or four concave 
spheroids (2CS, 4CS). The curvature criterion is d 2 Q/d\ Aj 2 < in 
the shaded region, indicating an instability toward a nonuniform su- 
perfluid phase. 

does not exist, and it is possible to map the grand-canonical 
and canonical phase diagrams onto each other as shown in 
Figs.^a) andgb). 



We emphasize that even though all of these gapless su- 
perfluid phases have order parameters with the same s-wave 
symmetry, they can be classified with respect to the topol- 
ogy of their zeros. It is well-known in other contexts that 
such changes in the topology leave strong evidences on the 
thermodynamic quantities such as the atomic compressibil- 
ity, spin susceptibility and momentum distribution, signaling 
a topological quantum phase transition. For instance, in nodal 
p-wave superconductors and superfluids for which the order 
parameters are k dependent, while the excitation spectrum 
= lAkp has gapless points of zeros in k space 

(which are determined by the nodes of | Ak| = 0) when /i > 
on the BCS side, it is fully gapped when /i < on the BEC 
side, showing a Lifshitz-type phase transition in the BCS-BEC 
evolution which occurs precisely at fi = 112311241 . The pre- 
sented phase diagrams clearly demonstrate here that it is pos- 
sible to simulate a much richer variety of gapless superfluid 
phases and numerous possible ways of phase transitions be- 
tween them with active experiments. 

Conclusions. To conclude, we analyzed the effects of both 
in- and out-of-plane Zeeman fields on the BCS-BEC evolution 
of a Fermi gas with an ERD coupling fields at zero tempera- 
ture. Depending on the combination of these gauge fields, 
we found novel gapless superfluid phases that can be distin- 
guished with respect to the topology of their zeros, e.g. one 
or two doubly-degenerate spheres, two or four spheres, two 
or four concave spheroids, or one or two doubly-degenerate 
circles. We noted that these phases are distinct because, even 
though their order parameters all have the same s-wave sym- 
metry, changes in the topology leave strong evidences on 
some of the thermodynamic quantities, signaling a topolog- 
ical quantum phase transition. 
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